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• ' We study a recently discovered family of 1/8-BPS supersymmetric Wilson 

■ loops in = 4 super Yang- Mills theory and their string theory duals. The oper- 
. ators are defined for arbitrary contours on a two-sphere in space-time, and they 

^ were conjectured to be captured perturbatively by 2d bosonic Yang-Mills theory. 

In the AdS dual, they are described by pseudo-holomorphic string surfaces living 
^ ' on a certain submanifold of AdS^ x S^. We show that the regularized area of 

' these string surfaces is invariant under area preserving diffeomorphisms of the 

boundary loop, in agreement with the conjecture. Further, we find a connection 
between the pseudo-holomorphicity equations and an auxiliary cj-model on S"^, 
which may help to construct new 1/8-BPS string solutions. We also show that 
the conjectured relation to 2d Yang- Mills implies that a connected correlator of 
two Wilson loops is computed by a Hermitian Gaussian two-matrix model. On 
the AdS dual side, we argue that the connected correlator is described by two 
disconnected disks interacting through the exchange of supergravity modes, and 
we show that this agrees with the strong coupling planar limit of the two-matrix 
model. 
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1 Introduction 

Exact results for non-protected operators in supersymmetric field theories which are 
valid for all values of the coupling constant A are typically beyond the reach of present 
techniques. For Af = 4 super Yang-Mills with SU{N) gauge gipup; AdS/CFT duality 
[l|-[3|, combined with the insights coming from integrability ^-[qI, allows to obtain a 
wealth of non-trivial information on the spectrum of the theory. Nonetheless, exact 
results which interpolate between weak and strong coupling regimes are rare. 
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A notable exception is the case of the half BPS circular Wilson loop 10l-ll2l|. For this 
non-local operator an exact answer, valid for all values of A and A^, is known. Indeed, 
it turns out that its expectation value is fully captured by a simple Hermitian matrix 
model with Gaussian potential. This remarkable result has been rigorously proved. 



using localization techniques, in 13| . It is clearly important to generalize this example 
and find a richer class of observables which preserve less supersymmetry but still allow 
for an exact answer as this will provide new interpolating functions between weak and 
strong coupling regimes. In 0, Il5| such a generalization was proposed by constructing 
a large new family of supersymmetric Wilson loop operators in the gauge theory. These 
operators are defined for arbitrary contours on a 5*^ subspace of Euclidean spacetime 
Mf^, and couple to three of the six scalars of A/" = 4 super Yang-Mills in a suitable 
way which depends on the spacetime contour. Generically these loops preserve only 
two supercharges (more precisely, two combinations of Poincare and superconformal 
supercharges). The half BPS circular loop is recovered as a special example when the 
contour is a great circle of S^. 

Of particular interest is the restriction to the sub-sector of loops which lie on a 
great S'^ inside 5''^. The restriction enhances the preserved supersymmetry to four 
supercharges. The dynamics of this sub-sector was conjectured to be fully captured 
by a reduced non-supersymmetric model, i.e. by bosonic 2d Yang-Mills on S"^, with 
a sp ecific prescription which amounts to discarding the contribution of 2d instantons 



15l . Il6| . According to this conjecture, the expectation value of the loops in A/" = 4 
super Yang-Mills can be computed by studying the analogous observables in 2d Yang- 
Mills, and restricting to the "zero-instanton sector". As shown in 174191], this sector 
is again governed by the Hermitian Gaussian matrix model, with a simple rescaling 
of the coupling constant by a factor which depends on the area of the loop on 5*^. In 
particular, for a loop along a great circle, one recovers the well known matrix model for 
the half BPS circle. A first indication in support of the conjecture was given at weak 
coupling in [l5|, [l6|, by observing that the combined "gauge field-l- scalar" propagator 
appearing in ladder diagrams effectively reduces to the propagator of pure 2d Yang- 
Mills in a suitable gauge with the WML prescription for the pole regularization. Further 



However, in 21 



non-trivial perturbative evidence was presented at 2-loops in |2d . |21 
a disagreement was found at order for the connected correlator of two Wilson loops 
along two latitudes of 5*^. 

In principle, it should be possible to prove (or disprove) the conjecture using lo- 



calization of the A/" = 4 SYM path integral, in the same spirit of [13 
step toward this goal has been taken in 



22 



A substantial 



23 , where it is shown that, due to the 



four supercharges preserved by the loops, the path integral indeed localizes on a 2d 
Higgs- Yang-Mills theory, which can be argued to be perturbatively equivalent to 2d 
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Yang- Mills. However the computation of the associated one-loop determinant in the lo- 



calization procedure appears to be technically challenging and was left open in [22| [23 
In the case of the half BPS circular loop, the determinant was in fact shown to be triv- 
ial 131]. If the same holds for the more general loops on S'^, the conjecture would be 



essentially proved. On the other hand, if the one-loop determinant is not trivial, it will 
deform the 2d theory with additional interactions which are not present in pure bosonic 
2d Yang-Mills. All results we present in this paper fit naturally with the original form 
of the conjecture, hence trivial one-loop determinant, but they do not strictly speaking 
exclude the other possibility. A non-trivial one-loop determinant may explain the dis- 



crepancy found in [21[ for the connected correlator. However, we point out that such a 
deformation of the 2d theory, if present, is quite constrained by the available data: it 
should not contribute up to 2-loops in perturbation theory and should be irrelevant in 
the strong coupling limit for arbitrary loops (due to the results presented in this paper 



and in [l5|, ll6|), and moreover it should not contribute to all orders for a loop at the 



equator, which is the half BPS circle for which the conjecture is rigorously proved 13 



Another possible explanation for the disagreement found in [21| is that there could 
be subtle issues in the choice of regularization scheme, which is not guaranteed to be 
compatible with the superconformal supercharges preserved by the loops. 

At strong coupling, Wilson loops in the fundamental representatioii^ correspond to 
macroscopic strings in AdS^ x 5*^ ending on the loops at the holographic boundary 



32, l33j. Their expectation value is obtained by computing the string action, after 
suitable regularization. From a purely mathematical point of view, AdS/CFT duality 
therefore instructs us to find the minimal area worldsheet with the correct boundary 
conditions and to calculate the associated area. For a loop of arbitrary shape, this 
computation is prohibitive. On the other hand, the conjectured relation with 2d Yang- 
Mills suggests, if correct, a very simple answer which is valid for a wide class of Wilson 
loops. A preliminary confirmation of the 2d Yang-Mills conjecture at strong coupling 



was given m 



15l . Il6l | by studying the string configurations dual to a loop at a constant 
latitude angle and to a loop comprising two arcs of longitude. 

One of the key-features of 2d Yang-Mills is the existence of a large symmetry group 
corresponding to the area preserving diffeomorphisms. To substantiate the 2d YM 
conjecture at strong coupling, we will prove that the regularized string action is indeed 
invariant under the action of this large symmetry group. To this end we will evaluate 
the variation of the on-shell string action as a function of the deformation of the 
boundary loop. In the calculation we will take advantage of the fact that the string 



^Wilson loops in higher rank representations can also be studied by replacing fundamental strings 
with D-branes with worldvolume electric flux, or more generally with the corresponding backreacted 



geometries, see e.g. |24426| |.|27 h31 |. 
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worldsheets have a simple geometrical interpretation as surfaces calibrated by an almost 
complex structure J of AdS^ x S'^, and live on a AdS^ x 5*^ sub manifold of the latter 6d 
space 151]. The boundary of AdS'^ coincides with the two-sphere on which the gauge 
theory loops reside. 

Supersymmetry usually drastically simplifies the understanding of string a-models 
by reducing the dynamical problem to the search of solutions of a set of first order 
equations. For the class of loops here considered, the first order equations are the 
"pseudo-holomorphic" equations associated to J. Despite this great simplification, 
finding explicit solutions to these equations remains rather challenging. It is therefore 
of great help, as we show in this paper, that the solutions of the pseudo-holomorphic 
equations are automatically solutions of a much simpler model, i.e. of a non-linear a- 
model on an auxiliary S^. If suitable boundary conditions are satisfied also the converse 
is true so that, given a solution of the auxiliary a-model, the full string configuration in 
AdS'i X S"^ can be reconstructed in a rather simple way. This approach offers therefore 
a new and insightful way to characterize the supersymmetric solutions of the string 
(T-model. For instance, by performing a reduction a la Pohlmeyer [s^], known solutions 
can be reinterpreted as solitons of complex sinh-Gordon. We will be also able to prove 
that any "stable" solution to the pseudo-holomorphic equations is accompanied by an 
"unstable" one with equal and opposite regularized area, generalizing the analogous 
result for the case of the latitude loop considered in 35||. This is in nice agreement 



with the strong coupling asymptotics of the Hermitian matrix model, where one finds 
two saddle points with equal and opposite exponents. 

The auxiliary a-model allows, in principle, to construct new solutions with greater 
ease. Indeed, as an application, we will find the string configuration which is dual to 
an infinitesimal deformation of the latitude solution. This "wavy" -latitude solution 
can be obtained by studying the linearized perturbations around the latitude solution 
and provides another check of the 2d Yang-Mills at strong coupling. We also find a 
simple generalization of the latitude solution, which has the topology of a cylinder 
and physically corresponds to a subdominant contribution at strong coupling to the 
connected correlator of two coincident latitude loops with opposite orientations. 

We then move on to study the case of connected correlators of two distinct (non- 
intersecting) Wilson loops on S"^. We show that the structure of the pseudo-holomorphic 
equations implies that supersymmetric solutions with topology of a smooth cylinder 
connecting two distinct latitude loops do not exist. Building on this result, and by 
using a "deformation" argument, we will argue that the same conclusion holds for two 
non-coincident generic loops on S*^. This result implies that at strong coupling the con- 
nected correlator is described by a "singular" cylinder consisting of two disconnected 
disks which interact by exchange of supergravity modes. 
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In the second part of the paper, we test our string theory conclusions against the 
Hermitian Gaussian matrix model which follows from the conjectural relation to 2d 
Yang-Mills. In particular, we show that the restriction to the zero-instanton sector 
implies that correlators of arbitrary 1/8 BPS Wilson loops on S"^ should be computed 
by an Hermitian Gaussian multi-matrix model. In general, exact results in 2d Yang- 



Mills with gauge group U{N) can be written as unitary multi-matrix integrals [36N46 
As noted in [12], discarding the contributions of the 2d unstable instantons amounts 
to neglecting the global structure of U{N), i.e. to approximate the gauge group with 
its Lie algebra, namely the (anti)-hermitian matrices. For the case of the two-loop 
correlator, we explicitly solve at large the corresponding two-matrix model and 
obtain the strong coupling asymptotic expansion. The saddle points appearing at 
large A precisely agree with the bulk expectation of disconnected disk solutions. Of 
course, the two-matrix model predicts a concrete set of subleading corrections to the 
exponential saddle points. It would be important to check those corrections from the 
bulk, by computing the precise amplitude for the supergravity mode exchange (as well 
as quantum fluctuations of the disks themselves), especially in view of the disagreement 
found in l2lll for the connected correlator. 



The paper is organized as follows. In Section 2 we will provide a brief review of 
the pseudo-holomorphic equations which govern the dynamics of the loops at strong 
coupling. In Section 3 we will prove the invariance of the string action under area 
preserving diffeomorphisms at the boundary. In Section 4 we will derive the auxiliary 
(T-model on and discuss the associated Pohlmeyer reduction. In Section 5, we will 
construct the "wavy" -latitude solution. In Section 6 we will focus on the coincident 
latitude correlator, while correlators between generic non-coincident loops will be ana- 
lyzed in Section 7. In Section 8 we discuss the matrix model computation for correlators 
of coincident loops, derive the Hermitian Gaussian multi-matrix model for correlators 
of non-coincident loops and solve explicitly the case of the two-matrix model. 



2 Review of pseudo-holomorphicity equations 

The Wilson loops introduced in are defined on a subspace of Euclidean 

spacetime M.^ which we can take to be given by x'^x'^ = 1@. They couple to three of the 
scalars of the Af = 4 gauge multiplet as 

W{C) =tTPexp ^ (^A^dx^" -i^^'^ , /i = 1, . . . , 4, i = 1, 2, 3, (2.1) 

^By conformal invariance, we can take the boundary sphere to have radius 1. It is straightforward 
to restore the radius dependence if desired. 
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where (t)j = 2a^^^jX^dx^ are the left-invariant one-forms on S'^ (we follow ll5i in the 



conventions for the SU (2) invariant one-forms j^l. The coupling to the scalars is such 
that the loops are 1/16-BPS supersymmetric operators, preserving two combinations 
of Poincare and superconformal supercharges. More details on the definition of the 
loops and on our conventions in the gauge theory will be given in Section 8. Here we 
will content ourselves with reviewing their dual string description. 

Since the loops couple to the three of the scalars, it is natural to consider an 
AdS^ X 5"^ subspace of AdS^ x . We write the metric on this subspace in the Poincare 
patch as follows 

ds^ = ^dx^'dx^ + z^dy'dy' , 

(2.2) 

z-^ ^ VY ■ 

This is related to the standard metric in the direct product form by changing variables 
as 0* = z?/*, where 6'* is a unit vector parameterizing S'^ . 

It was shown in [15^] that the string worldsheets dual to the general Wilson loops 
on S'^ reside on the AdSi,^ x S*^ subspace defined by 

xi'x^ + z'^ = 1. (2.3) 

Note that the boundary of the AdS^ is the on which the Wilson loop operators 



live. Furthermore it was shown in pj| that the dual string worldsheets are pseudo- 
holomorphic surfaces with respect to an almost complex structure J defined on this 6d 
subspace, i.e. the string embedding X^(r, a) satisfy the equations 

n^.X^ = ^ge^.d^X'' , (2.4) 

where M = 1, . . . , 7 and X^'^ = {x^,y^) are AdS^ x 5*^ embedding coordinates con- 
strained according to (12.31) . Here ^ is the determinant of the worldsheet metric, and 
our conventions for the antisymmetric e-symbol are such that era = ~1 (we denote by r 
the coordinate parameterizing the loop, and a the coordinate normal to the boundary). 
The explicit components of the almost complex structure J read 

K = z'al,y\ J':, = z'a;X = -z'J\. J'j = -zh^.^v' ■ (2.5) 

It is an interesting fact that the expression for J is closely related, up to factors of z, 
to the almost complex structure of which is inherited from octonion multiplication. 

It is possible to explicitly show that solutions to (12.41) automatically satisfy the 
equations of motion and the Virasoro constraints for the AdS^ x a-model. The 



•^An analogous construction in Minkowski spacetime, in which is replaced by the hyperbolic 



manifold H3, has been discussed in 47 1 
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reduction of the dynamics to first order equations is of course a consequence of su- 
persymmetry. In fact, by studying the K-symmetry condition for the string, it can be 
shown that the pseudo-holomorphicity conditions (12.41) are equivalent to the require- 
ment that the string solutions preserve the same set of supersymmetries of the dual 
gauge theory operators 

Given the almost complex structure J^, it is natural to use the metric to define 
a corresponding two-form Jmn- It then follows that pseudo-holomorphic worldsheets 
are cahbrated with respect to Jmn'- their area, or equivalently their on-shell action S", 
is given by 

S = ^[j. (2.6) 



27r 

Note that the two-form J is not closecQ and it is therefore not straightforward to obtain 
a general expression for the worldsheet area in terms of boundary data. However, it 



was shown in 15| that J admits a splitting 

J=Jo + dn. (2.7) 

where Q can be chosen in such a way that the exact piece dQ cancels the usual bound- 
ary term needed to regularize the infinite worldsheet area, namely the "Legendre trans- 



form" over the AdS radial coordinate z 49 



n= drzP,. (2. 

9S JdT. 



Therefore the regularized action Sreg = S — f^^ dr zP^ can be written as 



(2.9) 



In this paper we will be interested in the special subclass of 1/8-BPS Wilson loops 
living on a great S*^ C S^. The corresponding string solutions are obtained from the 
above general construction by restricting to the AdS^ x S"^ subspace given by x"^ = 0. It 
is then natural to use a notation for the coordinates such that x and y carry the same 
3d-vector index {x^,y^) — (x*,?/*). Accordingly, the components of the left-invariant 
one form simplify as o"^^ eijk- All pseudo-holomorphicity conditions in (12. 4p still 
hold, but the M = 4 component simplifies to 

9, (xV)=0, (2.10) 

which is readily solved by 

xy = c, (2.11) 

where C is an integration constant. 



*The non-closure of J points to a similarity with the case of a "generalized calibration" , see e.g. 
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3 Deformations of 1/8-BPS Wilson loops and in- 
variance under area preserving diffeomorphisms 



Consider an arbitrary curve C on S"^ corresponding to a 1/8-BPS Wilson loop operator. 
For any such curve, let us assume that we have found a corresponding dual string 
worldsheet X^'^(r, cr) which solves the pseudo-holomorphicity conditions fl2.4lfl Now 
consider an arbitrary small perturbation of the boundary loop C —>■ C + 6C such that 
the perturbed loop still sits inside S^. Correspondingly, the string solution in the bulk 
will be deformed to a new solution of the pseudo-holomorphicity equations 

X^\t, a) ^ X^'ir, a) + 5X^\t, a) . (3.1) 

For convenience, we will assume throughout that the solutions are in the conformal 
gauge ^/gg"^^ = 5"^, and to be concrete we can imagine that < r < 27r is a periodic 
coordinate parameterizing the loop C and < a < oo is the coordinate normal to the 
boundary, with the loop C sitting at a = 0. 

We would now like to determine how the on-shell string action varies as a function 
of the deformation of the boundary loop. In conformal gauge, the Polyakov action 
reads (for simphcity we set y/X/2n = 1) 

S = ^ f d-'ad^X^d^X'^GMN . (3.2) 

It is easy to see that the first order variation of the action resulting from (13.11) is a 
boundary term, since X^ is a solution of the equations of motion 

6S = ~ [ dT5X^^d^X^GMN. (3.3) 

Now we can use (12.41) to rewrite this as 

5S= [ drSX^^drX^JMN. (3.4) 

JdT. 

The integrand has to be evaluated in the limit o" ^ 0, so one needs to know the 
behaviour of the solution and of the perturbation close to a = 0. On general grounds 

^Though this is a reasonable assumption, we were not able to prove in general that solutions of 
(|2.4p exist for any given boundary loop C. An indirect argument based on holography is that, since in 
the gauge theory supersymmetric operators exist for arbitrary curves, we expect corresponding dual 
string solutions preserving the same supersymmetry to also exist. 
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one expects an expansion of the fornj^ 

x'{t, a) = xJo)(r) + aa;Ji)(r) + a^xl^^^r) + ... 

y'Ur) . . (3-5) 

Here x'^^^-^^t) parameterize the curve C on S^, and 2/(o)(''") will be related to the couplings 
6'*(r) to the scalar fields in the 1/8-BPS Wilson loop operator, see below. 

Plugging into (13. 4p the explicit components of Jmn (which can be obtained from 
(12. 5p by lowering an index with the metric) one can see that the Jij component does 
not contribute in the a ^ limit and the first order variation of the action can be 
written as 

dS = - [ dr {5x'dr{y^x^) + 5y'drX^x^) e^-fc . (3.6) 

By inspecting the expansion (13.51) one can see that only the first two non-trivial orders 
of X and 5X are needed in order to determine 5S. These can be obtained by solving 
perturbatively in a the pseudo-holomorphicity conditions and the constraints (12.31) and 
(12. lip . It is easy to see from the pseudo-holomorphicity conditions with M = 1,2,3 
that 

xl^ = 6x1^ = 0, (3.7) 
so the variation of the action is given by 



5S — — — J^^ dr (^5a;(o)9^(?/^Q^x^Q)) + 5?/(o)^(o)^0])) ^^J'^ 

dr (Sx\o)dr{yi.x'lo)) + 5?/[i)X^Q)xfo) ) eijk + 0{a) . 



(3.8) 



The leading two orders of y^ and Sy^ can be determined from the first two orders of the 
pseudo-holomorphicity conditions in the directions M = 5, 6, 7 and of the constraint 
()2.1ip . The leading order gives the condition 

1 • k 

i _ ^(0)^(0) .oqx 

which is nothing but the scalar coupling condition for the 1/8-BPS Wilson loop oper- 
ators [l6|], see the definition of the gauge theory operators in eq. (18. 5p (to obtain the 



above relation, we used the fact that |x(o)| = \d„z\„^Q, which generally holds for string 



^We assume that the AdSj, x S'^ coordinates (a;*(T, tr), z(t, cr), 0*(r, cr)) are analytic functions at 
cr = 0, i.e. they admit a Taylor expansion in positive powers of a. This is true for all known explicit 
solutions, though we do not have a first-principle prove of this assumption. The fact that y'(r, u) 
starts as — is just due to the definition = — and the fact that z vanishes at the boundary. 
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worldsheets dual to Wilson loops [49], and can be also seen to follow in the present 
case from the a-expansion of fl2.4p and (12. 3p ). Varying equation (13. 9p one then gets the 
corresponding expression for 5?/(o)- Plugging into the first line of (13. 8p one can see that 
it vanishes for arbitrary loop, so there is no divergent piece in 6S. Consistently with 
this fact, it is also possible to check that the boundary term which has to be added to 
the action in order to remove divergences, i.e. — J^j^zP^, does not give contributions 
at first order in 6X. Hence we have that 6S = SSreg (see also the discussion at the end 
of Section ij). 

At next order one gets the relations 



<^«ifc^(o)^(o)2/(i) = 

^(0)1/(1) = ^ 5 

where C is the integration constant in (12.111) . These are solved by 



(3.10) 



1/(1) = C'sjo) + a(r)i;Jo) , (3-11) 

where a(r) is an undetermined function, and correspondingly for the fluctuations 

= C6x\q^ + a(r)(5i;Jo) + 6Cx\q^ + 5a(r)i;Jo) . (3.12) 

Fortunately, the undetermined function a(r) as well as 6C and 6a turn out to be not 
needed in the evaluation of 6S. Plugging the above solutions into 6S and integrating 
by parts the first term in the second line of (13.81) . one finally obtains 

6S = —C / rfr^x^Q-ji'^Q-jX^Q-jejjfc . (3.13) 
t/ c 

Now one can recognize that the integral is precisely the first order variation of the area 
Ai of the loop C on S'^ (we define Ai to be the area of the region of S'^ bounded by 
the loop and including the north pole). One easy way to see this is for example to 
express the loop in polar coordinates {0{t), 4>{t)) where the expression for the area and 
its variation are 

Ai = dT(j){l - cos 6*) 

(3.14) 

8Ai = J dr sin e (^<i)6e -e6(j?j . 
Going to cartesian coordinates one can see that indeed 

SAi = y"rfr5a;(o)i(o)a;fo)eiifc, (3.15) 

so our final result for the variation of the action is 

6S = -C6A^. (3.16) 
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This proves that the string action is invariant under an arbitrary deformation of the 
boundary loop which preserves its area. Since the above deformation argument was 
carried out for arbitrary loop, it follows that any two loops with the same area will 
correspond to string worldsheets with the same value of the on-shell action. To fix the 
general functional form of the string action in terms of the loop area it is then sufficient 



to know its value for a simple class of loops, hke the 1/4-BPS latitudes of [35|[l5|, see 



eq. (I3.19P below, from which one infers (we reinsert here the explicit A dependence) 



Sreg = -^\/MA2, (3.17) 

where A2 = A — Ai {A = An is the total area of the sphere), in agreement with 
the conjecture that these loops are captured by the perturbative sector of 2d Yang- 
Mills theory on S"^, or equivalently by the Gaussian Hermitian matrix model with a 
rescaled 't Hooft coupling A — >• ^^^A, see Section [HI Note that in (13.171) we have 
written the result corresponding to the dominant contribution at strong coupling. It 



was pointed out in [35| in the case of the 1/4-BPS latitudes that for the same loop there 
are two supersymmetric solutions, one "stable" and one "unstable" , which have equal 
and opposite regularized areas: these correspond to whether the worldsheet wraps the 
smaller /larger region of the "dual" S"^ C with respect to the scalar coupling loop 
6'*(r). Intuitively one would expect this to be a general feature of arbitrary 1/8-BPS 
loops. In fact, in Section 142) we will give a simple argument that shows that for every 
given loop on the boundary S"^, there are two solutions to the pseudo-holomorphicity 
equations, which are related by C — > —C. From eq. (I3.16P one then sees that the two 
solutions should have equal and opposite regularized action. It is remarkable that both 
stable and unstable solutions appear as saddle points in the large A expansion of the 
Gaussian matrix model, as will be reviewed later. 

Note that, viewing (13.161) as a differential equation for S whose solution is (13.171) . 
the above calculation implies that the constant C does not depend on the shape of the 
loop and takes the value 

C = ±^%^, (3.18) 

where the +/— correspond to stable/unstable worldsheet as explained above. It would 
be nice to prove this relation by independent means. We have not been able to do this 
for the general case, but in the next subsection we will show that (13.181) can be proven 
to be correct in the limit of very small loops. Before doing that, we can show as a 
check that (13.181) agrees with the known solutions. For example, by looking at the ex- 



phcit latitude solution [35[[15[ (the choice of ± sign corresponds to the stable/unstable 
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solutions) 

1 tanhcTQCOsr ^ tanhaosinr o 1 

— 3^ — — 

cosho" ' cosho" ' cosh (To 

1 COST 2 sill''" 3 tanh((To ± 0") (3.19) 

^ z cosh(cro ± cr) ' ^ 2; cosh((To ± O") ' ^ z 

z = tanh (To tanh (T , 

and recalling that the parameter (Tq is related to the latitude angle by tanh (To = sin6'o, 
one gets indeed C = icot^^o? which agrees with fl3.18l) since Ai^2 = 27r (1 =1= cos^^o)- 
Another example of explicit solution is the 1/4-BPS string corresponding to a loop 



made of two half longitudes with an opening angle S 15|. In conformal gauge, the 
solution can be written a^ 

^ a sin aa sin a + cos aa cos a 2 ^ "^^^ ^'^ ~ ^'^ ^'^^ ^ 

cosh \J\ — a^T cosh v 1 — a^r 

= - tanh Vl - o?T , 7. = ^ '^^^ ^ (3.20) 



cosh \J\ — a^T 



^ sma(T 2 cosacT 3 



where the parameter a is related to the opening angle by a = and the range of the 
coordinates is —00 <r<oo, 0<(T<7r. From this solution we get C = x*?/* = yj=f , 
again in agreement with (13.181) after using Ai = 26 and A2 = 2(27r — 6). 

As a final remark, notice that the simplicity of the end result f l3.16p crucially de- 
pends on the structure of the pseudo-holomorphic equations in the special 1/8-BPS 
case. In the more general case of 1/16-BPS loops on S^, one could carry on a similar 
deformation argument with the result 



6S=- c?ry[i)(5xfo)X^o)<z. • (3-21) 



dT, 



Unfortunately in this case we cannot easily determine y^^^ as done in (13. lip for the 
1/8-BPS case, mainly because of the absence of the simple constraint x*?/* = C. 



3.1 The flat space limit 

Consider the limit in which the loop is very small compared to the radius of the sphere. 
Then the curvature of the sphere can be neglected and one ends up with a loop on flat 

^Here we write the "stable" solution. Analogously to the case of the latitude, there is a corre- 
sponding unstable counterpart, see Section ID 
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space. By studying the supersymmetry preserved by the corresponding gauge theory 
operators, it was shown in 15| IJ] that in this hmit one recovers the Wilson loops of 
Zarembo 50|] which preserve some fraction of Poincare supersymmetry. Correspond- 
ingly, we expect that eq. (12.41) should go in this limit to the pseudo-holomorphicity 
equations of 5l|], which describe Zarembo's loop in the bulk. Specializing to the case 
of loops on this limit can be systematically defined by setting 



X 



y 



■ex , 
e 



X = ex 
= 1,2,3, 



e z , 



(3.22) 



and taking the limit e — > while keeping the hatted coordinates fixed (note that x^ 
can be expressed in terms of and z using (12. 3p ). The pseudo-holomorphicity 

equations go in this limit to 



daX^ 



-2 o -2 

z eapOpy , 
0, 



daX^ = -feapdfjf , 



(3.23) 



which indeed are the same as the equations of 5l|] for the case of loops on (up 
to a trivial renaming of the y variables). The equation on the second line implies 
f = c, where c is related to C in eq. (12.111) by C = c/e. It is convenient to introduce 



coordinates on S'^ as 

e^ = z{f,-y\f) 
so that eq. (13.231) can be written as 

1 



dax' 
9' = 



(3.24) 



(3.25) 



cz 



where j„* = 6 daO^ — 6^da0 are components of the 5*0(3) current. This way of writing 
the equations displays a connection to an auxiliary S"^ a- model for the variables 6^, and 
this was used in [51| to show that the constant c only depends on the area of the loop. 
This can be seen in the following way. The area of the loop x on is given by 



2d 



d^o-eapdaX^djsx^ 



(3.26) 



Using an explicit parameterization 6^ = (sinw cos ly, sin f sin ly, cosf ) this becomes 



1 



TT 



= dwdv sin v cos v = — , 



(3.27) 



^The case of loops on should go exactly in the same way, and one will recover the equations of 
5l| for supersymmetric Wilson loops on R^. However the case of 1/16-BPS supersymmetric Wilson 



loops on cannot evidently be recovered as a limit of our loops on 
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where in the last equahty we have used that the worldsheet wraps half of the sphere 
(since z > 0). Since in the limit (13.221) the areas Ai, A2 defined before goes to 



Ai^e'A^d, A2^4n, (3.28) 
the formula (I3.27P implies that in this limit we have 

= ^ , (3.29) 

in agreement with the expected general relation (13.181) . 

The relation to the auxiliary S*^ a-model in (I3.25P is a special case of a more 
general relation to a S*^ a-model valid for the general Zarembo loops on ^l|. A 
similar connection to a a- model on an auxiliary holds for the case of our 1/8-BPS 
supersymmetric Wilson loops, as will be discussed in the next section. 



4 Pseudo-holomorphicity equations and cr-model 

In this section we show that, in the case of 1/8-BPS loops, there exist a correspondence 
between solutions to the pseudo-holomorphicity equations (12. 4p and solutions of an 
auxiliary cr-model. This may be useful to better understand general properties of 
these supersymmetric worldsheets and hopefully to find new solutions. 

For convenience, we will work again in the conformal gauge y/gg'^^ = 5"^. Using 
the explicit form of J in (12. 5p specialized to the case = 0, the full set of pseudo- 
holomorphic equations corresponding to the 1/8-BPS loops on can be conveniently 
written in vector notation as 

(4.1) 

Here the first two lines correspond to the M = /i component of (12. 4p and the last 
line to M = i. Note that these equations are not all independent. In fact, while (12. 4p 
are a set of fourteen equations, only at most seven can be independent, since = —1. 
So by this counting argument, one may expect that once the first two lines of (14. ip are 
solved, the remaining three equations should be automatically satisfied. This can be 
seen as follows. Define 

< = n^.^^ - ea,d,x^ 2) 



z da{xx y) = eapd/sx, 
x-y = C, 

—X X daX + z'^y X day = 
z^ 
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so that f ^ = gives the first two hnes of (14.11) and w^, = the last one. Using = — 1 
one can show that 

JWa = -Va - ^cpV^p , (4.3) 

SO if the first two hnes of (14. ip are satisfied one gets J^iWl^ = 0. Separating out the 
/i = 4 component and using the exphcit form of J this gives 

x^a^< = , = . (4.4) 

The first equation can have non-trivial solutions if the determinant of x^a^^j vanishes. If 
x'^ = 0, corresponding to 1/8-BPS loops, this in fact happens, and the zero-eigenvector 
is proportional to x\ So the first equation (14. 4p is solved by wl^ = CaX\ with Cq, some 
arbitrary function. But the second equation would then imply Cq, = 0, hence the only 
solution is = 0. 

This argument shows that it is sufficient to concentrate on the first two lines of (14. ip 
in order to solve the pseudo-holomorphicity conditions. The first equation implies 

d^{z''d^{yxx)) = 0. (4.5) 

Now define a 4-vector = (^, as 

^=zyxx, ^^ = VTTC^z. (4.6) 

Using the fact that 

{yx xf = y\' -(x-yf = \-l-C\ (4.7) 
where in the second equality we have used (12.30 and (12. lip , it is easy to show that 

ee = 1 , (4.8) 

so are coordinates on a round three-sphere. In terms of the eq. (14.50 becomes 

e^a^e - Cd^^^ = . (4.9) 

Multiplying this equation with and using C,^^^ = 1, it can be seen to be equivalent 
to 

d^^^ + da^''d^^''^^ = 0. (4.10) 
These are precisely the equations of motion of a a-model with action 

^ = 1 1 {d^i^d^i^ + A(e V - 1)) , (4.11) 

where we have introduced a Lagrange multiplier that enforces the condition ^^^^ = 1. 
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Note that since z > 0, the solution is constrained to he in the upper hemisphere of 
the S^. It is not difficult to work out the appropriate boundary conditions at a ^ 0. 
First of all z, and so has to vanish at cr = 0, and then by expanding the 3- vector ^"^ 
close to cr = in terms of x* and one concludes that the boundary conditions are 

eX=o = 0, eU=o = ^. (4.12) 

F{0)| 

So the appropriate solution is a worldsheet ending on the S"^ equator of this auxiliary 
S^, where it describes a loop parameterized by the unit tangent vector to the field 
theory loop x^'^^y 

So we have shown that to every solution to the pseudo-holomorphicity equations 
there corresponds a solution to an auxiliary a-model with boundary conditions 
fl4.12p . Clearly one would like to use this relation in the opposite direction and map 
the problem of solving the pseudo-holomorphicity conditions to the problem of find- 
ing classical solutions to the a-model. In principle, once a solution to (14.101) with 
the appropriate boundary conditions is given, the full solution x^,y^ may be be recon- 
structed as follows. First, the can be determined by solving the elementary first 
order differential equation, which follows from the first line of (14. ip 

d^x' = -^y^^^c^ {^%e - cd^e) . (4.13) 

As a remark, note that the right-hand side of this equation is proportional to a com- 
ponent of the conserved 50(4) current j^^ = i^d^i^ — i^d^i^- Once are known, 
the can be determined algebraically using the identity 

fxyxx=(l — z^)y — Cx , (4.14) 

which gives 

y = 0f X e + CSj , (4. 15) 

with z = vl + C^. The last line of (14.10 will then be automatically satisfied follow- 
ing the argument from (14.21) to (14.41) . 

In practice, however, it is not straightforward to apply successfully this algorithm 
because it is not easy to keep under control the boundary conditions on AdS^ x S"^ from 
the point of view of the auxiliary cx-model. In particular, once a solution on is 
given, it is not guaranteed that after integrating (I4.13P the resulting solution will satisfy 
the constraint x"^ + z'^ = 1. Besides this, there is also a self-consistency condition of the 
reduction to 5''^, namely = (see the definition eq. (I4.6P ). which is not guaranteed 
to hold in general from integrating (14.13p . It would be nice to understand if there is an 
efficient way to implement those constraints directly at the level of the variables. 
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As an example for which the algorithm can be efficiently applied, consider the 
following simple solution of the cr-model 

4 , i ( sinr COST \ 

e' = tanha, f= ^,_^,0 . 4.16 

\ cosho" cosher / 

This is a minimal surface inside S*^ which ends along a great circle of the equatorial 
S'^ . Geometrically, it is itself half of a great two-sphere inside S*^, as can be seen from 
the induced metric 

ds" = ^ (dr^ + da^) . (4.17) 

cosh a 

Solving (14.131) and (14.151) one easily finds that this solution generates the 1/4-BPS 
latitude in (13.191) . Note that, consistently with the fact that the unit tangent vector 
to a latitude is independent from the latitude angle, all 1/4-BPS latitudes correspond 
to the same solution on S^, the difference between them only comes from the choice of 
the integration constant C in (14.131) . This constant is fixed by imposing the constraint 
x^-l-z^ = 1, up to an overall sign choice. This sign ambiguity only affects the y*, see eq. 
(I4.15p . and corresponds to the stable and unstable solutions in (13.191) . The existence of 
stable and unstable worldsheets for each loop is a general fact which follows from the 
structure of the pseudo-holomorphicity equations, as will be explained in Section 14. 2[ 

As a remark, notice that not all solutions of the auxiliary a-model have an 
interpretation as minimal surfaces. This is only true if the satisfy the conformal 
gauge constraints on 

dredri"" - d^^^d^^^ = , dr^^d^^^ = . (4.18) 

In this case, ^^{t, a) is a minimal surface and the action (14. lip gives its area. But the 
Virasoro constraints for the string a-mo del on AdSs x S"^, namely T^p^+T^^ = 0, do not 
in general imply the vanishing of the stress tensor of the auxiliary a-model. In fact, 
it is possible to show from the pseudo-holomorphicity equations that the constraints 
(I4.18P are satisfied if and only if the Virasoro constraints on AdSs x S'^ are separately 
satisfied, T^^^ = = T^^. This happens for example in the case of the latitude solution, 
and hence the corresponding solution on has a minimal surface interpretation, but 
it is clearly not the most general case. 

An example of a relevant solution which does not have a geometric interpreta- 
tion as minimal surface, is the one corresponding to the two-longitude solution. The 
conformal gauge constraints on are not met in this case, since the two-longitude 
solution on AdSs x only satisfies T^p^ + T^i^ = 0. The explicit solution on is 
given by 

= — cos acr tanh Vl — o?t , = sin aa tanh Vl — a^r , 

^3 COSO" ^4 smo" (4.19) 

cosh a/1 — cAt cosh Vl — a?T 
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Solving fl4.13p and (14.151) one can recover the corresponding supersymmetric solution 
on AdSs X 5*^. The constant C is fixed by the constraint + = 1 to be C = ±-^^=5. 
The choice of positive sign gives the solution in ( 13.20p . while the choice of negative sign 
gives a different solution which we interpret as the unstable counterpart of (I3.20p . This 
unstable solution, which has not been previously discussed in the literature, is related 
to the stable one in (13.201) by ^ — ^f'^^^ f?- It is easy to see that this solution 
has the same boundary conditions as the stable one, but wraps the "complementary" 
region on S"^ C compared to (13.201) . It can be also checked by explicit calculation 
that its regularized action is equal and opposite to the one of the stable solution, as 
expected. 

To conclude this section, let us also mention that the reduction to the a-model 
allows to rewrite the general expression for the regularized action, eq. (12.91) . in a po- 
tentially more suggestive form. Since z oc it follows from the 5''^ a-model equations 
of motion that 

Sreg = j^'^O = ^ j d'^T [dj'dj' - d^^C^d^C^) . (4.20) 

This equation is for example useful as it allows to easily repeat the deformation argu- 
ment of Section [2] directly at the level of the regularized action. Since the equations 
of motion for the 9^ G S"^ and the G are separately satisfied, it follows that the 
variation of Sreg is a boundary term 

SSreg = "J / dr {59%9^ - 5^^^^^^) . (4.21) 

Going through the same steps as in Section [3l one finally finds that SSreg = —C6Ai, 
in agreement with the fact stated earlier that the regularizing boundary term does not 
affect the deformation calculation of Section [31 



4.1 Pohlmeyer reduction and Complex-sinh-Gordon model 



It is well known j34 , that given a a-model on we can perform a reduction a la 

Pohlmeyer which reduces the dynamical variables leaving only the "physical" degrees 
of freedom. This kind, of reduction was also recently considered in various string related 
contexts, e.g. in 



55-64 



In this section we will apply the reduction to the auxiliary 
sigma model on the S^. This will shed a new light on some of the solutions of the 
pseudo-holomorphic equations. Hopefully this might also help in gaining some insight 
on the general properties of the solutions and to construct new ones. 

Exploiting the classical conformal invariance of the a-model we can always bring 
the components T^z and Tzz of the stress tensor into the form 



dz^^dz^^ = d-z^d-ze = 



dz = d^- idr ,dz = d^ + idr 



(4.22) 
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where is a constant. To obtain the "reduced" theory we now introduce the fields 
and 9 through the following non-linear field redefinition 

cosh 20 = (9,^^95^^, 



-f^-%e sinh^ = -dieK"" , (4.23) 
fx^d.e sinh' = ^dl^^K^ , 
where the vector K is defined as follows 

K'' = eABCD^''d,fd,e- (4.24) 

The use of hyperbolic functions, compared to the trigonometric ones as in the usual 
Pohlmeyer reduction, is convenient in this Euclidean worldsheet setting for keeping 
the variable real. Via this procedure we have effectively eliminated one of the fields 
of the cr- mo del and reduced the theory to a two-dimensional one with 6 and as 
dynamical variables. If we introduce the complex variable 

^ = sinh0 exp(-ze/2) , (4.25) 

the equations of motion for the original cr-model are satisfied provided that ip satisfies 
the complex-sinh- Gordon (CshG) equation 



^^-^^^ - '^*Y+W ^ + 1^1) = • (4-26) 

When 6 is constant this equation reduces to the sinh-Gordon equation for the variable 



<9f9^0 + sinh(20) = . (4.27) 

After Pohlmeyer reduction, the longitude solutioij^ (14.191) satisfies the CshG equa- 
tion with n = a and can be written as 



y/l — a"^ exp(— zo") 

Remarkably, ipiong is formally analogous to the soliton solution of complex-sine-Gordon 
theory (CsG). Indeed the solitonic solution of the latter theory can be written as 
followfl 0, E 

exp(m sin at) , , 

i^soi = cosa—^ (4.29) 

cosh(/i cos ax) 



^We note that the latitude solution ()4.16p is singular from the point of view of Pohlmeyer reduction 
since in this case /i = 0. 
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The most general soliton solution can be obtained by "boosting" eq. (|4.29p . 
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where a is a constant parameter. The dyonic magnon 5^ (which is a generahzation 
of the giant magnon solution 66|) is a sohton of CsG which corresponds to ipsoi- The 
longitude solution can be recast in a form equivalent to (14.291) by taking fi = a and 
cos a = In the case of ijjiong the role of "time" t is taken by the variable 

a while r is the analogue of the spatial coordinate x. 

The CshG theory can be thought as a generalization of sinh-Gordon with an addi- 
tional internal U(l) symmetry. The CshG solitons are charged under this symmetry. 
The longitude solution ipiong carries a U(l) charge which can be computed as follows 

Q = z f d^tMjI^t^ = 4arctanh(v/r^). (4.30) 

Remembering the relation a = (tt — 5)/7r between a and the longitude angle 5, we see 
that the charge of ipiong grows monotonically with 5 and vanishes for 5 — > 0. Up to a 
rotation, the latitude solution corresponds to the singular limit of infinite U(l) charge. 



4.2 Unstable solutions 

In the examples of the latitude and two-longitudes loops we have explicitly seen that 
there are two different solutions, corresponding to the same loop at the boundary, which 
preserve the same supersymmetry but have equal and opposite regularized action. In 
this section we argue that this is a general property of all 1/8-BPS loops on S"^. 

Suppose we have a solution {x\y^) to the BPS equations (14.11) . corresponding to 
a given supersymmetric loop on S"^. Then we can construct a new solution (x*,y*) by 
keeping the same while changing the as 

f = y'- 2C-^ . (4.31) 

Note that the AdS radial coordinate z = ^/y~^ is left invariant by this redefinition of 
the I/*. Since xxy = xxy, it is clear the first line in (14. ip is still satisfied. Further, 
we have that = —C, so the second line of (14. ip is also satisfied, with the new 
integration constant simply given by C = —C. As argued at the beginning of Section 
m the last line of (14. ip follows from the first two, so we have shown that (x*,^*) is a 
new solution of the pseudo-holomorphicity equations. Moreover, the new solution has 
the same boundary conditions as the original one, since the are the same, while the 
second term in (I4.3ip is subdominant close to the boundary, so that zy^\„=o = 2;y*|o-=o- 
Since C —C, we expect from eq. (13.160 that the solutions {x\y^) and {x\y^) should 
have equal and opposite regularized areas. It would be interesting to also prove this 

^^This choice corresponds to a — tt/2 purely imaginary. 
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directly by comparing the explicit expressions (12.91) for the two solutions. While the 
interpretation of this pair of solutions as stable and unstable worldsheets is natural, 
it would be also interesting to directly study the quadratic fluctuations around these 
saddle points, and show that for one of the solutions there are tachyonic modes in the 
three directions orthogonal to the S"^ C 5*^ (generalizing the observations of 35| for 
the 1/4-BPS circular loop). 



5 Linearized perturbations around circular solutions 

In this section we obtain an approximate 1/8-BPS solution by solving the pseudo- 
holomorphicity equations for linearized perturbations around the latitude solution. 
The reduction to the a-model obtained in Section H] turns out to be quite useful for 
this purpose. The results of this section give an explicit confirmation to the general 
argument of Section [3], and they also provide an example of how to use the cr-model 
to generate AdS^ x S'^ solutions. 

It is convenient to take the following conformally fiat metric on 5'^ 



^ ^ V + . (5.1) 



This form of the metric is related to the embedding coordinates by 
. ^ f 2pcos0 2psin0 2C 1 - - e \ 

^ Vi + p' + c^'i + p' + c^'i + p' + c^'i + p' + cv ■ 

The (T-model action in these coordinates reads 

S = - f d^CT 1 ^"^^"^ + P^(9«0(9a0 + de^CdaC 

27 (i+p^ + ef ■ 

We now define the fiuctuation fields around a solution (p, 0, () as 

p = p + ^(i+p' + C')<^i 

= 0+^-(l + p2^C')</'2 (5.4) 
z p 

As it is often convenient, we have rescaled the fluctuation fields by the inverse vielbein 
so that ifi carry flat indices. We have also introduced a formal parameter e to keep 
track of the perturbative expansion. 
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To work out the linearized equations of motion for the fluctuations, one should now 
expand the action at quadratic order and read off the corresponding kinetic operators. 
Let us specialize to the case in which the solution to expand around is the one given 
in eq. (14.161) . corresponding to the 1/4-BPS latitude 

-p = e-% = r + |, C = 0. (5.5) 

Dropping total derivatives which do not contribute to the equations of motion, the 
action at quadratic order then turns out to be 

S{2) = Y y '^^^ ^(•^ctV'i)^ + (<9ctV^2)^ + {dr'^i - tanh 0-^22)^ + (<9tV52 + tanhcrcy^i)^ 

1 2 

72— + 'A) + daVsdaVa 71—^3 

cosh a cosh a 

(5.6) 

This can be viewed as the action for a free scalar (the transverse fluctuation ip^) and 
a free vector field (the two longitudinal fluctuations) living on the S"^ hemisphere with 
metric (I4.17p . To see this, note that the two longitudinal fluctuations can be combined 
into a vector field A"" = {(pi, (f2), where a is a 2d tangent space flat index. The covariant 
derivative with respect to the metric (I4.17P acts as 

VaA'' = do.A- + u^\A\ (5.7) 

where the only non-zero component of the spin connection is uj^'^^ = — tanhcx. Then 
the quadratic action can be written covariantly as 

^(2) = y y d'a^ {g-^V^A^V^A^ - A'^A^ + g^^d^cp.df^ifs - 2ifl) , (5.8) 
or, integrating by parts and dropping total derivatives 

^(2) = - Y / ^'^v^ (^"^^' + + V + 2)ips) . (5.9) 

The equations of motion for the fluctuations 

V^A- = -A\ VV3 = -2<^3, (5.10) 

are now recognized to be the equations satisfied by vector and scalar spherical harmon- 
ics on 5*^ 

VVj'" = -[/(/ + 1) - 1] rj'" , VV''" = -/(/ + 1)F'™ , (5.11) 
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with Z = 1 in both cases. The solutions are well known and can be written in terms of 
Legendre polynomials 

r) = P,'"(tanhcr)e^'"" , yi'"(cr, r) = e«^9^y""(cr, r) , (5.12) 

with |m| < Specializing to the case / = 1 and taking real linear combinations, one 
gets the following solutions of (IS.lOp 

cos T sin T 

(^3 = aotanhcr + Oi h &i — , (5.13) 

cosh a cosh a 

and 

(fi = = cosh ada(p , (p2 = A'^ = — cosh adr'^ , 



COST ~ smr 
if = do tanh a + ai — h bi - 



(5.14) 



cosh cr cosh a 
where ao, ai, 6i, ao, fii , &i are arbitrary constants. 

To obtain the perturbed solution in terms of the coordinates x\ on AdS^ x S"^, it 
is first convenient to go back to the embedding coordinates for which the solution 
to linear order in the perturbation is 

.1 sinr , . , , 

^ = e ( sm T tanh cry^i + cos T(p2) , 

cosh a 



= h e (cos r tanh a^pi — sin rv?2 

cosher 

= tanh cr 



COST 



cosher ' (5.15) 



5 

'4 _ ^ V'l 



cosh cr 

From here one sees that the boundary condition ^^ocz— >Oatcr = fixes fio = 0, 
while all other constants are still arbitrary. To get the solution for the x*, one has then 
to integrate 

d.x' = ^,A_ 6„^ {^%e - cd^e) . (5.16) 

Since we are considering a perturbation of the latitude solution, the integration constant 
C appearing here can be written in the formic 

C = cot 00 + £7 + C(e^) , (5.17) 



^^There is a small subtlety here since we are looking for solutions on the hemisphere with 
< cr < cxD rather than the full sphere — cxd < cr < oo. The restriction on the quantum number m 
comes from requiring regularity on the whole sphere. If we only need regularity at ct — s- cx) and cr — > 0, 
then additional solutions of the form P™ (tanh cr)e*™'^ with — oo < m < —I would be allowed (here 
P™(tanh cr) has to be thought in general as an associated Legendre function and it is not a polynomial 
for m < —I). However one can see a posteriori that including these extra modes does not lead to 
a solution satisfying the correct boundary conditions and the constraint a:^ + = 1, so we do not 
include them here. 

""^•^ Again we choose a specific sign of C corresponding to the stable latitude solution. One can 
analogously obtain a perturbed unstable latitude solution. 
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for some constant 7. One can now solve the above differential equation to linear order 
in e and fix the resulting integration constants so that the constraint x'^ + z"^ = 1 and 
the consistency condition = are satisfied. The latter constraint fixes Oq = 0, 
and one obtains the perturbed solution for the function of the free parameters 

7, ai, bi, di, bi and the latitude angle Oq. The explicit expression is a bit lengthy and 
we do not write it out here. The loop at the boundary x^q-^ = x*(r, 0) is some kind of 
"wavy" latitude whose area can be computed to be 

A,= rf /(o)^(o) ^(0)^(0) =2n{l- cosOo - e^sm^Oo) + 0{^) . (5.18) 

Jo 1 + ^(0) 

Note that the area only depends on the constant C consistently with the general relation 
f l3.18p . though the shape of the loop and of the worldsheet clearly depends on all other 
constants. 

Having determined the and the x*, the are explicitly obtained as explained 
above from eq. (14.151) . The resulting solution (x*, y*) satisfies the pseudo-holomorphicity 
equations up to quadratic order in e. 

The regularized area of the solution can be computed from fl2.9p and the result to 
linear order in e turns out to be 

Sreg = sm Oq {1 - e^f sm Oq COS Oq) + Oie"^) . (5.19) 

Comparing with (15.181) . one finds agreement with the general formula Sreg = — V^i^2- 



6 SUSY solutions for coincident latitudes 

The reduction of the holomorphicity conditions to the a-model can be used to obtain 
a new super symmetric solution, which we interpret as corresponding to the correlator of 
two coincident latitudes with opposite orientations. This can be obtained by adapting 



to our case a solution on given in [51 



Start with the following parametrization of 

= {cost] cos (j), COST] sin (f), sin rj cos -ip, sin T] sin 'ifj) , (6.1) 
and take the circular symmetric ansatz (the shift by 7r/2 in is for later convenience) 

TT 

(f) = T + -, n = r]{a), ij = ij{a). (6.2) 

The (T-model action reduced on this ansatz reads 

S d'^a + sin^ ri^p'^ + cos^ rj) . (6.3) 
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There are two conserved quantities 

= sin^ r]tp' 

1 (6.4) 
L = - [r] + sm rjip — cos rj) . 

We now restrict to the case in which ip' = and take ip = 7r/2. Then 77 (cr) is determined 
by 

rj''^ = — sin^ r] , (6.5) 

where we have defined = 2E + 1 > 0. If e = 1 this differential equation is elementary 
and it is solved by sin 77 = tanhcr, which is the solution (14.161) . and thus corresponds to 
the single latitude solution. For general e the above differential equation can be solved 
in terms of elliptic functions 

^ = ± //"^ . 2 = • (6.6) 

Jo ^/e^~sm^r| e V ^ J 

The behavior of the solution is different depending on whether e > 1 or e < 1 and the 
undetermined ± sign must be chosen accordingly. 

For e > 1, the derivative 1]' never vanishes, and //(a) increases monotonically 
(corresponding to the positive sign in (16. 6p ) from 77(0) = to 77(0"/) = vr, where 
(7/ = ^F^^n,-^). Correspondingly sin 77 (and hence the AdS coordinate z), starts 
at zero at a = 0, reaches a maximum sin 7/ = 1 at a = ^ and then vanishes again at 
a = af. The topology of the solution is thus that of a cylinder, since there are two 
boundaries. 

For e < 1 the derivative t]' vanishes at rj = t/max = arcsine, therefore one chooses 
the positive sign in (16. 6p for < a < ^ and the negative one for ^ < a < aj, where 
now the final value of a is given by (Xf = (arcsine, The coordinate r] therefore 
starts at zero, reaches a maximum 77max < f and then turns back to reach again r] = 0. 
Also in this case sin 77, and hence z, has two zeroes, so the topology is again that of a 
cylinder. 

Through equations (I4.13P and (14.151) we can now reconstruct the full supersymmet- 
ric solution in AdS^ x 5*^. Firstly it is clear that 



sm T] 



(6.7) 



Integrating (I4.13P one gets 



1 , , v^e^ - sin^ T] 2 / ^ V ^2 - sin^ 7] . 3 a/1 + - 

X = sie) ■ cos r , x = sie) , — smr, x — 



VTTc^ ViTc^ ViTc^ 

(6. 
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where s(e) = 1 for e > 1 and s{e) = sign(^ — cr) for e < 1. For all values of e, the curve 
at both boundaries of the cylinder is the same latitude circle with angle related to 

C= ^''-f^r (6.9) 
sm Oq 

Note that for e < 1 the solution exist only as long as e > sin^o- 
The coordinates can be obtained from 

y= ^\ i^^^^s X f + cs] , (6.10) 

+ cos"^ 7] \ sm?7 / 

where ^ = {S} ■, ■, S,^) ■ The explicit expression is a bit lengthy and we do not write 
it out here. It is possible to explicitly verify that the full solution (x*,y*) correctly 
satisfies all the pseudo-holomorphicity equations (14.11) . In both the e < 1 and e > 1 
cases it can be seen that the loop ends on opposite latitude circles on the S"^ C 5*^, 
namely on zy^ = sin^^o at a = and on zy^ = — sin^^o at a = cxf. This difference in 
the scalar coupling is consistent with the supersymmetric boundary conditions and is 
due to the opposite orientations of the two coincident latitudes in space-time. 

Though in both cases the solution describes a correlator between coincident lati- 
tudes, the behavior of the solution in the bulk is qualitatively different for e > 1 and 
e < 1. In particular, for e < 1 the (x^,x^) circle shrinks to zero size in the bulk while 
the {y^,y^) never does, and the opposite happens in the e > 1 case. In the limiting 
case e = 1 both circles shrink to zero size at the midpoint. 

It is not difficult to evaluate the regularized area of this solutions by using eq. 
(12. 9p . Plugging in the solution and evaluating the integral one finds that, both for 
e > 1 and e < 1, the contribution from the region < o" < ^ is exactly cancelled by 
the contribution coming from the other half ^ < c < af. So the conclusion is that in 
both cases 

Sreg = 0. (6.11) 

Note that the value of the action, as well as the boundary conditions, do not depend on 
the parameter e. We interpret this parameter as a modulus which should be integrated 
over with the appropriate measure in the string partition function, but we do not 
investigate this further here. 

Besides the connected solution described above, one can also consider "degener- 
ate" solutions for correlators of coincident latitudes using the known 1/4-BPS latitude 
(I3.19p . A solution with the topology of the cylinder can be obtained by joining two 



"'^^For simplicity we pick here the positive sign of C. There is an analogous solution corresponding 
to the choice of negative C. 
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disks of the type (13.191) . one for each of the two coincident loops, by a thin zero-area 
tube. Physically, this should be interpreted as the exchange of light supergravity modes 
between the two worldsheets (note that the two worldsheets are not coincident because 
the scalar coupling is different due to the opposite orientation). The area of the so- 
lution is just going to be equal to the sum of the two areas, while the amplitude for 
the supergravity mode exchange will contribute as a prefactor of the exponential in 
the string partition function. By coupling this way the stable/unstable versions of the 
solution (13.191) . one obtains worldsheets with regularized area 

Sreg = — 2V^sin6'o , Sreg = 2V^sin6'o , Sreg = . (6.12) 

The solution obtained by gluing the stable and unstable worldsheet can be thought as 
the limiting case e — 1 of the above family of connected solutions with zero regularized 
area. The first choice with Sreg = — 2-\/Asin6'o is clearly the dominant contribution at 
strong coupling. However we will see that all three possibilities appear as saddle points 
at strong coupling in a matrix model computation. 



6.1 The coincident correlator as soliton of sin(h)-Gordon 

Also for the solution (I6.2l) - (l6.6p contributing to the coincident latitude correlator we 
can apply the Pohlmeyer reduction and interpret it as a solution of CshG. In the present 
case it reads 

cn(ecr, e~^) , „ , s. , ^ 

V'cor. = Ul— ^ exp i-teo/2) , (6.13) 
Ve — 1 

where cn(ecr, e~^) denotes a Jacobi elliptic function. Since ^ = is constant, this solu- 
tion has zero U(l) Noether charge. Given a chargeless solution ip = sinh0 exp(— i^o/2) 
of CshG, (I) is automatically a solution of sinh-Gordon and in the case of the correlator 
we havqlfl 

= arcsinh (6.14) 

which is indeed a solution of eq. (14.271) with /i^ = — 1. This expression is formally 
analogous, upon analytic continuation, to a "kink-wave" periodic solution of the sine- 



Gordon equation (see e.g. |56l. l67l|) which can be written, in the static limit, as 

arcsin(cn(fca;, — )) . (6.15) 

The connection with sine-Gordon solitons can be seen directly, even without performing 
Pohlmeyer reduction, by noticing that the variable t] appearing in the circular sym- 
metric ansatz satisfies the sine-Gordon equation d'^rj + sinrj cosr] = 0. The solutions 



15 



For e < 1 it is convenient to add i7r/2 to (|6.14p to keep real. 
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relevant for the correlator are obtained by inverting (16.61) and are of the forno 

7] = am(ecr, 1/e^) , (6.16) 



which is another form of the kink-wave solution present in the literature, see e.g. [& 
Using rj as variable allows to safely perform the single latitude limit e ^ 1 to obtain 

7] = arcsin(tanh(cr)) (6-17) 

which is the familiar expression for the kink solution of sine-Gordon. 

7 Correlators of non-coincident loops 

It is natural to ask whether there are supersymmetric string solutions describing con- 
nected correlators between two generic loops. In this section we will provide evidence 
that a connected smooth configuration does not exist, the only possible surface being 
a degenerate one, obtained by joining the worldsheets of the individual loops by an 
infinitesimal tube corresponding to an exchange of supergravity mode^. We will first 
start by showing that a connected worldsheet does not exist for the case of the corre- 
lator of two latitudes loops. Using a "deformation" argument we will then argue that 
supersymmetric connected solutions do not exist for loops of arbitrary shape. 

Let us begin by considering the connected correlator of two 1 / 4-BPS latitudes with 
latitude angles, say, 9i and 62. To find the associated string we should look for a 
solution to the pseudo-holomorphicity equations with the topology of a cylinder whose 
boundaries are given by the two latitude circles. Because of circular symmetry, we can 



take a periodic ansatz on AdSs x S of the type considered in [69 

X* = tanhp (sin u COS'?/', sinw sin-?/', cosm) , z 
= cosh p (sin v cos 0, sin v sin 0, cos v) , 



coshp (7.1) 



witflisl 



p = p(cr) , u = u{(t) , V = t>(cr) 
-0 = r + a(o") , (j) = T + b{a) . 



(7.2) 



^^For e < 1 one has to invert the sign of 77 for cr//2 < a < <jf. 

^^In the case of two 1/4-BPS latitudes, the same conclusion was previously reached independently 
by Nadav Drukker 0. 

^^One may in general take arbitrary winding numbers — k\T + a{a) , — + (5{a), but the 
supersymmetric boundary conditions require k2 — ki, and we take the single wrapping case ki = 1. 
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The relevant solution should have x^{a) = tanhpcosu varying continuously from 
x^(0) = cos^i to x^lcxf) = cos ^2- However it is not difficult to see that the pseudo- 
holomorphicity equations do not allow such a solution. If we take the first line of eq. 
(14. ip with a = T and plug in the above ansatz with circular symmetry we obtain 

d^x^ = . (7.3) 

So x^ has to be constant and we fall back to the single latitude solution (or the coin- 
cident latitudes of the previous section). We then conclude that there is no connected 
super symmetric solution describing the correlator of non-coincident latitudes. The only 
relevant supersymmetric worldsheets will then be degenerate cylinders constructed by 
joining two standard latitude solutions (13.191) by a zero area tube, as described at 
the end of Section [6l corresponding to the propagator for the exchange of supergrav- 
ity modes between the two worldsheets. To leading order in the a' expansion, the 
connected correlator will be given by the exponential of the regularized area of this 
degenerate surface, which is just the sum of the areas of the two disks. Considering 
as before all possibilities of stable/unstable pairings, we then find four different saddle 
points with regularized areas 

Sreg = VX{± slu Oi ± slu 62) . (7.4) 

We will show in Section 18.31 that a Gaussian two-matrix model calculation precisely 
exhibits in the strong coupling limit the presence of these four saddle points. 

We now discuss the general case of a correlator between two generic distinct loops 
of areas Ai and A2. Let us suppose that a connected surface with cylinder topology 
joining these two loops does exist. We can then act with an area preserving deformation 
on the loops to bring ourselves to the case where we have two latitudes whose angles 
6i are determined by Ai and A2. As a consequence, the initial connected surface will 
be deformed into the one describing a latitude-latitude correlator. As proved before, 
this final surface must correspond to a disconnected worldsheet. This therefore implies 
that, while deforming the loops, we have moved from an initial connected worldsheet 
to a disconnected one. This fact signals the presence of a discontinuity, i.e. a phase 
transition 70|-[72|, at some intermediate configuration as we interpolate from the initial 
loops to the latitude ones. This phase transition, which is typically of first order, should 
manifest itself as a singularity of the free energy of the string. To leading order in the 
supergravity approximation, the string free energy is given by the worldsheet area. 
But by assumption we kept the area of the loops constant during the deformation 
and, for our class of loops, by the arguments of Section [3] this in turn implies that 
the worldsheet area of the string representing the correlator is also constant as we 
deform the loops. As a consequence the free energy cannot change and therefore no 
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singularities can possibly appear during the deformation. This leads to conclude that 
our initial assumption was inconsistent and that as a consequence the connected surface 
joining two loops with the topology of two latitudes does not exist. Hence, we conclude 
that the connected correlator is described at strong coupling by the four saddle points 
fl7.4p . with sin^j ^ ^^'^^2"^ (the dominant contribution being, of course, the one 
corresponding to the choice of two negative signs). 



8 Matrix model computation 

In this section we study connected correlators of Wilson loops in the zero instanton 
sector of the two-dimensional Yang-Mills on a two-sphere. Two-dimensional Yang-Mills 
theory is a well-known and explicitly solvable theory 361-146 . 

The exact partition function of the 2d Yang-Mills on arbitrary Riemann surface 
S can be obtained by gluing the basic building blocks of three types: the partition 
function of 2d Yang-Mills on a disk, on a cylinder, and on a sphere with three holes. 
Each of these partition functions depends on the boundary conditions of the gauge 
field, and in the gluing process one integrates over the boundary conditions. For a 
boundary with topology of a circle 5"^ the boundary conditions for the gauge field 
modulo gauge transformations are parameterized by the holonomy of the gauge field 
around S^. The holonomy takes values in the gauge group G. Hence, if S is glued 
from the basic building blocks along K circles we get the partition function of the 2d 
Yang-Mills on S represented by an integral over K copies of the group manifold G. 
Of course, one can insert Wilson loop observable on any k-th gluing circle, which will 
be simply represented by inserting a character of the corresponding k-th integration 
variable in the integral over K copies of G. We will refer to the partition function of 
2d Yang-Mills represented in this form as the unitary A'-matrix model. 

In the simplest case of a sphere with possible insertion of one non-self-intersecting 
Wilson loop we get a unitary 1-matrix model. A Lie algebra q of Lie group G can 
be viewed as a perturbative approximation to G. For G = U {N) the Lie algebra g is 
represented by Hermitian matrices multiplied by V 

The connection between the zero-instanton sector of two-dimensional Yang-Mills on 



S"^ and the Hermitian matrix model was shown in ITHlOl. |73| . Intuitively, discarding 
instantons in 2d Yang-Mills exactly corresponds to discarding the non-trivial topology 
of G, and replacing G by its tangent space at the unity, i.e. by the Lie algebra of G. 
The exact partition function of 2d Yang-Mills on a disk, which is a heat kernel for the 
standard Laplacian on G, is replaced by the heat kernel of the standard Laplacian on g. 
But the heat kernel on g is a simple Gaussian function. Therefore one relates the zero- 
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instanton sector of 2d Yang-Mills on a sphere with one non-self-intersecting Wilson 
loop to the Hermitian matrix model. For concentric Wilson loops on a sphere the 
argument goes in the same way, and we explicitly derive the corresponding Hermitian 
multi-matrix model in Section 18.31 



8.1 Conventions 

We take Euclidean space-time and assume that the gauge theory action is normal- 
ized as 

Sym = j d^xti Qf^.F^'^ + D^'I'.D^*^ + ^[<l>„ $,][$\ <^^] + (8.1) 

where the symbol 'tr' is a negatively defined invariant bilinear form on g. For a simple 
group G the choice of the bilinear form is unique up to a scaling. For the gauge group 
G = U{N) we take 'tr' to be the trace in the fundamental representation, hence the 
notation. Assuming that {T°}, a = 1 . . . dimg is the orthogonal basis on q, and A"" are 
the real coordinates on g in this basis, i.e. A = A"'Ta, and assuming the normalization 
tr T"T'' = —^6ab the propagator in the Feynman gauge is 

(^»K(y)) = If ^73^- (8-2) 

Our conventions are such that the gauge field A takes value in the Lie algebra of the 
gauge theory, e.g. in the anti- Hermitian N x N matrices if the gauge group is U{N). 
The covariant derivative is D = d + A and the gauge field strength is the two-form 
F = dA + AAA. 

For a G-gauge bundle over the space-time M and two points xi,X2 E M connected 
by a path C, the holonomy of a gauge connection A is a parallel transport of the gauge 
bundle fiber from the point X2 to the point xil^ 

G(xi,X2,C) = Pexp / A^dx^. (8.3) 

A holonomy G{xi^X2,C) takes values in the gauge group G. For a closed contour C 
and a representation R of the gauge group the Wilson loop operator is the /^-character 
of the holonomy around C 



lyfi(C) =tr^Pexp lA^dx^. 



Jc 



•4) 



^^The path ordering is defined in such a way that if is any point on C then G{xi,X2) — 

G{xi, Xm)G{Xm, X2)- 
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In A/" = 4 SYM it is natural to generalize the definition of Wilson loop operator 
by adding couplings to the scalar fields of the A/" = 4 multiplet. This allows to define 
operators which globally preserve some fraction of the superconformal symmetries. In 
particular, here we consider the 1/8-BPS Wilson loops for contours C restricted to a 
unit sphere S*^ C M^, which are obtained as a special case of the general loops in fl2.ip 
by restricting to a great two-sphere inside S'^ . They couple to the three scalar fields as 

Wr{C) = tiR Pexp ^{Aj + ieijk'^'x'')dx^ . (8.5) 

These Wilson loops can be thought of as holonomies of the complexified connection 
Ac = A+i*^ (here * is a two-dimensional Hodge star on 5*^) in the 2d Higgs- Yang-Mills 



theory [22, Um 



lyfi(C) =tr«Pexp jfi,. (8.6) 



In |14|-|16[|. it was conjectured that the expectation value of a non-self-intersecting 
Wilson loop (18. 5p can be computed in the Gaussian matrix model with coupling con- 
stant which depends on the area on 5*^ enclosed by the loop 



{Wr{C)) = ^ [ [d<p]e ^^''^'tr^e^ (8.7) 
^ -he 



where the partition function Z is, of course, fixed by saying that {Wji{C)) = 1 if i? is 
the trivial representation. In ( 18. 7p the integral is taken over the imaginary slice of the 
complexified Lie algebra Qc, i.e. over Hermitian matrices in the U{N) case. The Qmm 
is the area-dependent effective coupling constant in the matrix model 

By A = 47rr^ we denote the total area of 5*^, and Ai, A2 are the areas inside and outside 
the loop C on 5*^; they satisfy Ai + A2 = A. For example, specializing to the case of 



the latitude at the polar angle 6 measured from the North pole we have 28|, |35l 

~g'yMiO)=9YM^^'0, (8.9) 

SO 9 = 71 /2 corresponds to the equator and gives the standard 5-BPS Wilson loop 



so 9 




10, 


11 



To streamline notations for the matrix model computation we denote 

2 N 



9mm 



— ; X = gi,j^jN = Aiy, (8.10) 
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then the matrix model has the form 

{Wr{C)) = i / [d<P]e~^^'^'^\iRe^ (8.11) 
^ J Is 

so V has a meaning of effective 't Hooft couphng constant in the matrix model. 

The arguments which lead to the matrix model conjecture are exactly the same 
if we consider several Wilson loop on top of each other in arbitrary representations, 
hence we have 

{Wr,{C)Wr,{C)...Wr,{C)) = ^ f [d<P]e-^"'^" tTR.e^ tin.e* . . .tiR.e'^ . (8.12) 

^ J Is 



8.2 The connected correlator of two coincident Wilson loops 

We are going to compute the connected correlator of two Wilson loops on the same 
base contour C 

{Wr.Wr,),,^^ = {Wr.Wr,) - {Wr,){Wr,) . (8.13) 

We specialize to the case where i?i = / is the fundamental representation ofG = U (N), 
and i?2 = / is the dual (complex-conjugate) representation to -Ri. This corresponds to 
two coincident loops of opposite orientation^. The large N expansion of the two-point 
correlator in the Hermitian matrix model has been previously considered in 7^. For 



results at finite A^, see 75 



The A;-point resolvent function is defined as 

GiZi, ...,Zk)= iV^-'(tr . . . tr ^-)eonn . (8.14) 

zi- (p Zk- (p 

This function is a generating function for all connected correlators of the form 

(tr0*Hr0^^..tr0*'=),onn, (8.15) 

since 

oo oo 

Giz,, ...,Zk)= N"-^ ^r'"'' • • • (tr tr . . tr 0*'=),o„„ . (8. 16) 

Another generating function W{xi, . . . ,Xk) is 

W{xi, ...,Xk) = A^'-2(tr e^i^ . . . tr e^^<^)co„n = 



oo oo 

5Z ■ ■ ■ 5Z ^1' • • • 4' (-^ -7 tr <P'' tr . . tr 0^^)eo„n 



(8.17) 



^"The case of same orientation corresponds to R2 — f and it is also easily obtained from the 
following calculation. 
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The function W{xi, . . . ,Xk) can be obtained from G{zi, . . . , Zk) by the inverse 
Laplace transform in all variables zi, . . . ,Zk 

W{x,, ...,Xk) = jdz^...j dzkG{zu ;2,)e-i-i+--'=-^ , (8.18) 

where the integration contours encircle the point z = oo for all integration variables. 

The fc-point resolvent function can be expanded in inverse powers of N, which is 
the same as the topological expansion over the genus h of Riemann surfaces associated 
with Feynman ribbon graph diagrams 

oo 

G{z,, ...,Zk) = J2 Ghizi, Zk)N-''' . (8.19) 

h=0 

The leading order term h = corresponds to the planar diagrams. 

T he g enus zero fc-point resolvent can be computed recursively from the loop equa- 



tions 



76l478l|. The 1-point and 2-point genus zero correlation function are 



Go{zi, Z2) 



Z1Z2 - 4z/ 



(8.20) 



2(^1-^2)2 I ^zj - Au^zl - 4z/ 



To compute Wilson loop correlators we need to perform the inverse Laplace trans- 
form. For the 1-point function the result is 

where Ji is the modified Bessel function of the first kind. For a single Wilson loop in 



the fundamental representation we take x = 1 and recover the well-known result 10 



The large u asymptotic can be extracted from the integral representation for the 
modified Bessel function /n(a;)0 

/„(x) = ^ [ e"™^(^) cos(r2^) dO . (8.23) 



The modified Bessel function J„ is related to the Bessel function J„ of imaginary argument 
/„(x) = (-i)"J„(ix). 
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At large x we can use the semi-classical approximation. There are two critical points 
6 = and 6 = tt. Expanding in 6 near 6 = and taking the integral perturbatively we 
get 



.X = . e 



1 ^ ^^ 1-4^2 1 



(I6n^ - 40n2 + 9) + ...^ . (8.24) 



V 8x 128x2 

Expanding in 6 near 6 = tc we get get 

Notice the fl8.25l) is imaginary for real x because we make the expansion around the 
unstable critical point@. 

In this way we get the large u expansion for a single Wilson loop 
Woix = = e'^ (^^7^ - 320F^5/4 " i024v/?z.V4 + " " " ) 

V 2Viz/3/4 ~ 320Fz/5/4 + i024v^i/7/4 + 

It is remarkable, as first pointed out in jssi, that the two saddle points corresponding 

to e^"^^ = have a bulk interpretation as stable and unstable supersymmetric 

worldsheets with equal and opposite regularized discussed in Section [^^21 

Now we consider the two-point connected correlator. It is not trivial to make inverse 
Laplace transform of the function Go{zi,Z2) in f l8.20p because it does not factorize. 
However, one can make a trick differentiating Gq{zi, Z2) over the parameter v and then 
integrating back [7^. Indeed, 

a.Golzi, Z2) = ^,2 _ Auy/^ ^^-^^^ 

we see that dvGo{zi^ Z2) as a sum of two factorisable terms. Making inverse Laplace 
transforms of the factors in (18.271) 

1 xli{2^x) 



.2 _ 4j,)3/2 2^ 

z 



[z -^ur' -v<^ /8_28) 



[z^ - Avfl^ 



xIo{2^/ux) 



22 



The full asymptotic series is [79[ p. 223 



^ (-l)fc(n,fc) e-"^("+^)" ^ {n,k) 
"(-)-^l. (2,). + V2^ 1.(2.)- 

where {n,k) := {—l)''{n + l/2)k{—n + l/2)fc/fc! and Ofc is the shifted factorial: ao = l,afc = a{a + 
1) . . . (a + fc — 1) for an integer fc > 0. The ± sign must be taken + for — 7r/2 < arga; < 37r/2 and — 
for — 37r/2 < argx < tt/2. 
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we get 



(8.29) 

After integration over v with boundary condition W^i^z^^ Z2)\v=q = we get 
/ X _ ^^1^2 {h (2v^xi) Ji {2^x2) + /q (27^x2) h (2x/^xi)) 

KI/ol^Xi, X2j — ; (^O.dUj 

Xi + X2 

To evaluate the correlator of coincident Wilson loops with opposite orientation we need 
to take the limit x\ = —X2 = 1 of the above expressioij^. We get 

Woixu X2)\,,=-„ = -x^v (/o (2xv^) 2 + I2 (2xv^) h {2x^) - 2h {2x^) ^) . 

(8.31) 

The large u asymptotic expansion gives 

l^o(xi,X2)U,=-.,=i,™ = e^v^ ("Y^^ - - 20487rz/3/2 + 

^""''^^ (leV^ " 128^ ^ 20487rz/3/2 + ' ' ' ' ^^'^^^ 
^ \ n IOttv^ ~ 10247rz/3/2 ^ ' ' 

Notice that there are three contribution to the series corresponding to the exponents 
of 4:^/u, —Ay/u and 0. Indeed, as explained in Section [6] in the case in which the loop is 
a latitude on S"^, these contributions can be observed in the dual type IIB string theory 
in the AdS^ x background. They correspond to string solutions with regularized 
areas equal to Sreg = ±2\/X and Sreg = 0, see eq. fl6.11l) - fl6.12l) . 

We also quote the small u expansion of the connected correlator for reference 

Wo{xuX2)U=.,,=i,u-,o = -y- \v' - - - ^/ + • • • (8.33) 

8.3 The connected correlator of two distinct Wilson loops 

Ar gum ents based on the localization to the two-dimensional theory 15, 16| and 13, 
231] lead us to the following i^-matrix model to compute correlators of K non-self- 



22 



23 



The choice of same orientations corresponds to = X2 = 1. ExpUcit expressions for this case 



can be found e.g. in [74 1 
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intersecting Wilson loops which have topology of K latitudes on S'^ 

{Wr,{C^) . . . Wr,{Ck)) = ^ / [dct>i] . . . [# 

^ J la 



_ K\ 

Is 

^^ii^^i+^i^2-M^+-+^^^(^k+i-M^+-+^{<t>Kf) (8.34) 



where is the coupling constant of the 4d A/" = 4 SYM, Ai is the area enclosed 
by the first loop, ^.+1 are the areas between the k-th and the k + 1-th Wilson, and 
Ak is the area enclosed by the last loop. There are K + 1-terms in this matrix model 
corresponding to the K + 1 effective propagators on the cylinders between Wilson 
loops. The first and the last propagator is contracted with the Dirac delta like wave- 
function corresponding to the the vacuum states on the North and the South poles of 
the sphere!^ 

To be more concrete, recall that the two-dimensional Yang- Mills on a cylinder with 
the action 

S2d = -i / d^x^tiF^f.F'^^ (8.37) 

in Hamiltonian formalism is equivalent to the quantum mechanics on the group mani- 
fold G with the Hamiltoniaio 

- H = ^A = ^C2, (8.38) 
4 4' ^ ' 



where A is the Laplacian on G, the C2 is the second Casimir operatoCj, and the time 
is measured in the units of area thanks to the usual area-preserving diffeomorphism 



invariance of the 2d Yang-Mills, see e.g. 36|-|46 



2" A "quick and dirty" way to get the matrix model (jOi)) is to recall that [40|, |4l|, |80| 2d YM 
on Riemann surface S localizes on classical configurations: cIa * F = 0. Also, in two dimensions, 
Yang-Mills theory essentially reduces to the abelian theory 40|,|4l|, so in our case we can consider 
to be a piecewise constant function on E with jumps supported on Wilson loops. The contribution to 
the 2d YM action coming from 'E,k,k+i (the region between fc-th and (fc + l)-th Wilson loop) is 

Sk^k+i = I Sx^iv{*Ff = ~^Ak,k+itrXl^+^ , (8.35) 

where Xk^k+i denotes ^F on Y.k.k+i- Using Stokes theorem on "Ek^k+i 



Ak,k+iXk,k+i ^ d^x^'^F^ F^ A- A 

'I^k,k + l ^Ck+\ ^Ck 



Ofc+1 



.36) 



we get (|8.34p and the relation = — ^^xal tjjis derivation we assumed that after the abelianization 



the gauge bundle is trivial, and this precisely means dropping out the unstable instantons [41 1 



4C 



^^In Euclidean space our convention is ^ = —H^ where \1/ is the wave-function. 
^^In our conventions the basis Ta on g is normalized as tr^;' TkTb = — \5ab and we choose the metric 
on g in this basis to be ^5ab- The Laplacian for this metric (Casimir operator) is then A = C2 = 



37 



The configurational space here is the holonomy of the gauge connection around 
a spacial circle, i.e. it coincides with the gauge group G. The wave-function if) is 
a function on G, i.e. we write ip{U) for U E G. The heat kernel K{U,l\t) for the 
Laplacian f l8.38p defines the wave-function of the theory at time t starting from the 



Dirac-delta function supported at the group unity at time t = [361 . 
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K{U,l-t) = J2dRexp{-^C2iR))xRiU), (8.39) 



R 



where xr is the character of U G G in representation R and the sum is taken over all 
irreducible representations of G@ More generally the propagator between two adjoint 
invariant states represented by Ui and U2 is 

KiU^,U2;t) = J2xRmxRiU2)expi~^C2iR)). (8.40) 

R 

The expectation value of K Wilson loops on a cylinder located at times -|- 
A12, . . . and with trivial initial and the final state is then 



[dU,] . . . [dUK] K{1, f/i; A) K{U^, U^, As) • • • K{Uk, 1; Ak) x (8.41) 

tr/ji Ui tr/j2 U2 ■ . . tiRi^ Uk ■ 

The integral in fl8.4ip is taken over K copies of the group manifold G. 

However, we can relate this integral to the integral over the Lie algebra of G up to 



the exponentially small instanton corrections [l7j. We consider G = U{N) gauge group 
for brevity, but the same arguments are applicable to any simple compact Lie group. 
Let 9i be coordinates on the Cartan subalgebra f) of g, and let diag(e''^^, . . . , e*^'^) be 
the image of (^1, . . . , 6]\f) G f) under the exponential map g ^ G. Clearly the integral 



^Cae" — 2^^eaea, where Ca are basis elements in g corresponding to T^. The constant C2{R) is 
defined as Ta{R)Ta{R) = ~C2{R)1. Hence A^-^ = -2C2(i?)*i? if is a matrix element of G in 
representation R. 

^^We remind that the characters for a compact simple group G are orthonormal 



1 



G 



volG 

and that the delta-function at the group unity is expressed as 

R 

where the sum is taken over all irreducible representations of G. 
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in fl8.4ip is invariant under adjoint transformations of Uk for k = 1 . . . K, and hence 
we can reduce it to the integral over diagonal unitary matrices. Let for z = 1 . . . 
be coordinates on the maximal torus T G G with T being the image of f) under the 
exponential map, and let tij for i ^ j he local coordinates on G in a neighborhood of 
T. Let Tij,i 7^ j be generators of g/f), so that tijTij G g/f) parametrize the normal 
directions to T. One can take A = diag(e*^i, . . . , e*^^), then W = e*^^"^^^ and 

U = WAW-^ . (8.42) 

In these coordinates the metric on G in the neighboorhood of T has the familiar form 



N 

i=l l<if^j<N 



^ ft -ft 

J2de^+ Y: (2sin(^))^<.. (8.43) 



Using this metric and the standard formula for the Laplacian 

Ag = -^d.i^g^'^d,) , (8.44) 

we can find that on the adjoint invariant functions on G represented as functions on T 
and denoted as K{9i, . . . , 6' at) the Laplacian acts as 



where 



J(^^) = n(2sin(^^)). (8.46) 



i<j 



Now notice that 



It is a very nice algebraic identity which can be proved using Weyl formula that the 
last term is a constant 

J^^W = -^,NiN-' - 1). (8.48) 

hence we get 

^gK = {JK) + ^N{N^ - l)K. (8.49) 

The addition of a constant to the Hamiltonian only shifts the zero energy level and 
does not change the correlation functions, so for the purpose of computing correlation 
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functions of Wilson loop operators we can drop it. Then we end up with a simple 
formula 

AgK = j-^ (JK) . (8.50) 

Repeating the same arguments for the Lie algebra q instead of G we can get Lapla- 
cian on the adjoint invariant functions on g 

= i-JoK) , (8.51) 

where 

MO) = ll{Oj-e^). (8.52) 

Looking at the Laplacian fl8.49p on G and the Laplacian (18.511) we notice that heat 
kernel equation for the functions Kq = JKq and = JqKq is the same usual heat 
kernel equation in the flat space 




(8.53) 



However, strictly speaking, the solution Kg on G and Kg on g is not the same. The 
difference is the global structure of the space. The exponential map g ^ G is a local 
isomorphism between g and G but, of course, not a global one. In particular, the points 
on the integer lattice (27ml, • • • , 27™^) in f) are mapped to the unity of the G. Hence, 
while Kg satisfies the heat kernel equation fl8.53p with a single delta-function source 
inserted at the point 9i = 0, the Kq satisfies the same heat kernel equation fl8.53p but 
with infinitely many images of the delta-function source inserted at the points on the 
lattice (27rni, . . . , 27mk),ni, . . . ,nk G Z. The contribution of the images for n 7^ is 
exponentially suppressed for g2d — ^ and corresponds to the instanton contribution to 
the 2d YM partition function. 

Rewriting the unitary matrix model f l8.4ip as integral over eigenvalues 9k for k = 
1 . . . we get (here Ui = UiiOi) ...Uk = Uk{Ok)) 

j [d^'e] . . . [d^'e] Kcil, 61; A,)J{ei)J{9i) Kg0i, 92-, A^2)J02f . . . KciOK, 1; Ak) x 

Xr{0i)---Xr{0k). 

(8.54) 

We argued that J{ei)KG0ij2,t)J{e2) is equal to Jo(^i)^g(^2, ^2, ^)^o(^2) up to 
the instanton corrections. In other words, the heat kernel Kq on invariant functions 
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on G multiplied the measure factors J{6) which appear in the reduction of the unitary 
matrix integrals to the eigenvalues, is equal, up to instanton corrections, to the heat 
kernel Kg on invariant functions on q similarly multiplied by the factors Jq which 
appear in the reduction of the Hermitian matrix integral to the eigenvalues. Hence the 
"zero-instanton" contribution to the correlator is 

XR{e^)...XR{eK). 



].55) 



Explicitly, for generic 6i and 62 we have 



M9,)K,{9,,92,t)M92) = const -i- J] (-l)'^exp (-^-M^ " ^(^"2))') , (8-56) 

where Sn is the permutation group of N elements. For a general Lie algebra the 
permutation group should be replaced by the Weyl group. It is clear that the right 
hand side of fl8.56p satisfies the usual heat equation on the flat space fl8.53p . Also, 
when we reduce the integral over the Lie algebra to its Cartan subalgebra there is still 
a discrete residue symmetry group left - the Weyl group. The summation over Weyl 
group in (18.561) projects to the Weyl invariant state. The heat kernel Kg{9i,92,t) is 
a positive function and is Weyl invariant (completely symmetric in 9u or 92i), while 
Jo{9i)Kg{9i,92,t)Jo{92) is Weyl skew-invariant. 

If the source is the origin of the Lie algebra, i.e. at 9i = 0, then the solution (I8.56P 
is degenerate and is not applicable. The correct solution is 

K,{l,92,t)M92) = const^Jo(^"'2)exp (^—1-(92Y^ ■ (8-57) 



For more details see e.g. [ITj and references therein. 

Replacing the integral (I8.55P over the eigenvalues of the Hermitian matrices by 
the integral over all Hermitian matrices we recover the Hermitian i^-matrix model 
chain (18.340 up to a factor of i = ^/—l in the definition of observables. In other 
words, the Wilson loops in the 2d YM with compact gauge group G are unitary, 
e.g. in the fundamental representation X/(^) = SjLi 6xp(27ri6'j) with 9j being real. 
The observables in the Hermitian matrix model which describe the supersymmetric 
Wilson loops of the c? = 4 A/" = 4 Yang-Mills are not unitary, but Hermitian, i.e. 
the fundamental Wilson loops are represented in the Hermitian matrix integral as 
X/(^) = XljLi ^^P(2^^j) "with the same convention that 9j are real. 
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One can equivalently formally represent d = 4 Af = 4 Yang-Mills SUSY Wilson 
loops by the usual unitary Wilson loop in the bosonic 2d Yang-Mills in the zero in- 
stanton sector and with the inverted sign of the 2d YM coupling constant 

9l = (8.58) 

[A = ATcr"^ is the total area of the 5*^). 

In the case of two Wilson loops we can write the two-matrix model action in the 
form 

Z2MM = j M0i][rf02]e-^*'"(^^('^^)+^^('^^)-^'^^'^^) , (8.59) 

with two quadratic potentials Vi(</)i) = V2(02)) = 1^202 the interaction 

coupling c. In terms of the original areas and the 't Hooft coupling constant A = gyM^ 
we have 

ai 







1 






h — 

A12 






1 




U2 


h -— 

A12 


A 


1 





.60) 



XA 



12 

We now solve in the planar limit the correlator (18.341) in the case of the two-matrix 
model and for choice of fundamental representation. Following 8l|] , define two resolvent 
functions 

1 

zi = ax + ai— (8.61) 

X 

Z2 = Oi — h a2X , (8.62) 

X 

where a,ai,Q!2 are coefficients determined from the algebraic equations 

^ (8.63) 

-V^{Z2) - Zi = X. 

c ca 
Solving these equations on a,ai,a2 we find 



a 



2 _ c 



aia2 — 



ai = —a (8.64) 
c 
ai 

a2 = — a . 

c 
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Now, inverting the equation (18.611) we get the resolvent function xi{zi) and inverting 
the equation (18.621) we get the resolvent function 0:2 (-22) 



xi = ^ [ zi± \/z'f - Aaai 
1 

2ao 



X2 = - — { Z2 ± \/ Z2 - Aaa-^^ 



(8.65) 



In 81j and 8^ it was shown that the planar limit of the two-point resolvent function 



for the matrices 02 



uj{zi^Z2) = {ii—^ti^—) (8.66) 



conn 



is given by 

To cast this formula in more symmetrical form with respect to the two loops we will 
denote ^ 

yi{zi) = xi{zi) 1/2(22) = — ^ (8.68) 

X2\Z2) 

SO 

uj{zi, Z2) = -9^1 log (1 - yi{zi)y2{z2)) , (8.69) 

where now 



^2(^2) = ^ ( ^2 + V - 4aa2 



^.70) 



To compute correlators of Wilson loops we need the loop generating function 

Vr(ti,t2) = (tre*i<^Ure*^'^^>^^^^. (8.71) 

The function ^^(^1,^2) is obtained from uj{zi, Z2) by the inverse Laplace transform 

W{h,t2) = (fdzidz2u;iz,,Z2)e'''^+'''\ (8.72) 

where the integration contours are sufficiently large to encircle the cuts in the zi and 
the Z2 plane. 

The choice of signs in (I8.70p and the integration contours is fixed by comparing the 
resolvent with the leading order Feynman diagrams computed directly from the matrix 
model fra]) . 
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8.3.1 Large A asymptotic 



So we get 



Now we will compute fl8.72p by the stationary phase method in the limit A ^ oo. From 
f l8.67p we have 

(1 - yi{zi)y2{z2)y 

It is convenient first to change integration variables in the integral fl8.72p from Zi 
to yi 

Zi = ayi + — i = l,2. (8.74) 
Vi 

(27r?)2 / (1 - yiy2Y 
We can integrate over the circular contours of radius > 1 

yi = r,e"*-\ r, = .p, 0, e [0, 27r) (8.76) 
V ot 

so 

^ ' ^ (27r)2 y ^2(1 _ rir2eWi+*<^2))2 ^ ^ 

Notice that for not coincident loops rir2 > 1 (see fl8.82p ) so we replace the integrand 
by the convergent series in inverse powers of (rir2) 

oo 

2/12/2 



Each term in the series factorizes into a simple product of Bessel functions. Using 
( K23\f we get 

oo 

W{h, t2) = 5^ r2(rir2)-"/„(2tiv/^)/n(2t2v^) • (8-79) 



n=l 



Let us denote 



(y.\(y.2 

p = rir2 = J ti = 2tiy/aai t2 = -2t2^Jaa2 ■ (8.80) 

We introduced the minus sign in the definition of t2 for convenience in order to specialize 
to the case of oppositely oriented Wilson loops for which t\ > 0,^2 < 0, so both ti,t2 
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are positive. Then 



n=l 



1 1 



1 



1 



Si v^2^ V(l + P) 



-P /e*i+*2_g-f2-fi^+, P 



1-p)^ 



g-tl+t2 ^ gil-t2 ) + . . . 



f8.811 



The resuh in terms of the areas Ai,A2,Ai2 and 't Hooft couphng constant A is 
easily obtained using (18.801) at ti = —t2 = 1 and (18.601) . (I8.64p . which give 



a 



XA1A2 
A^ 



ti = 2^aa.i = 2 



A(A-Ai)Ai 



A 



A 



t2 = 2^aa2 = 2 



X{A- A2)A2 



A 



A 



^2) 



P 



a\a.2 



1 + 



A 



12 



1 + 



A 



12 



A. 



Therefore, the prediction of the two-matrix model for the strong coupling behavior of 
the correlator is, at leading order. 



W{\, -1) exp 



A4Ai(A - ^) ^ , /A4A2(A - A2) 



A^ 



A^ 



(8.83) 



which matches with the bulk interpretation given in Section [7] of a singular cylinder 
made up of two disconnected disks joined by a supergraviton propagator. Moreover, 
the presence of four saddle points in (I8.8ip is also mirrored in the bulk by the four 
possible pairings of stable/unstable disk solutions. It would be very interesting to test 
the structure of the prefactors to the exponential saddle points in (18.811) by computing 
the supergraviton exchange on the string theory side. 
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